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ABSTRACT 

We present a new determination of the primordial helimxi mass fraction Yp, based on 93 spectra of 86 
low-metallicity extragalactic H ii regions, and taking into account the latest developments concerning 
systematic effects. These include collisional and fluorescent enhancements of He i recombination 
lines, underlying He i stellar absorption lines, collisional and fluorescent excitation of hydrogen lines 
and temperature and ionization structure of the H ii region. Using Monte Carlo methods to solve 
simultaneously for the above systematic effects, we find the best value to be Yp = 0.2565 ± 0.0010(stat.) 
± 0.0050(syst.). This value is higher at the 2a level than the value given by Standard Big Bang 
Nucleosynthesis (SBBN), implying deviations from it. The effective number of light neutrino species 
is equal to 3.68;^;?^; (2cr) and 3.80l:^;?g (2cr) for a neutron hfetime t„ equal to 885.4 ± 0.9 s and 
878.5 ± 0.8 s, respectively, i.e. it is larger than the experimental value of 2.993±0.011. 
Subject headings: galaxies: abundances — galaxies: irregular — galaxies: ISM — ISM: abundances 



1. INTRODUCTION 

The determination of the primordial "^He (hereafter 
He) abundance and of some other light elements such 
as D, •^He and ^Li, plays an important role in testing 
cosmological models. In the standard theory of big bang 
nucleosynthesis (SBBN), given the number of light neu- 
trino species, the abundances of these light elements de- 
pend only on one cosmological parameter, the baryon- 
to-photon number ratio 77. 

Because of the strong dependence of its abundance on 
r/, deuterium has become the baryometer of choice. The 
D/H measurements in damped Lya systems appear to 
converge to a mean primordial value log D/H = —4.56 
± 0.04, correspon ding to a baryon m ass fraction ri^/i^ 
= 0.0213 ± 0.001 (|Pettini et al.H2008t ). This estimate of 
ribh^ is in excellent agreement with th e value of 0.02273 
± 0.00062 obtained by Dunklcy et al.l ([2009) from anal- 
ysis of five years of observations with the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP). 

Although He is not a sensitive baryometer (Yp depends 
only logarithmically on the baryon density), its primor- 
dial abundance depends much more sensitively than that 
of D on the expansion rate of the Universe and on a pos- 
sible l epton asymmetry in the early universe (jSteigmanl 
I2OO6I I2007D . Thus, accurate measurements of the pri- 
mordial abundance of He are required to check the con- 
sistency of SBBN. 

However, to detect small deviations from SBBN and 
make cosmological inferences, Yp has to be determined 
to a level of accuracy of less than a few percent. The 
primordial abundance of He can, in principle, be derived 
accurately from observations of the He and H emission 
lines from low-metallicity H 11 regions. Several groups 



have used this technique to derive the primordial He 
mass fraction Yp, with somewha t differ ent results. In the 
most recent study llzotov et al.l ()2007[ ) based on a large 
sample of 93 spectra of 86 low-metallicity extragalac- 
tic H II regions (hereafter the HeBCD sample) derived 
Yp = 0.2472 ± 0.0012 and Y^ = 0. 2516 ± 0.00 1 1, us - 
ing lBeniamin et al.l (|1999ll2002[ ) and lPorter et al.l H2005f ) 
He I emissiviti e s, re spectively. On the other hand, 
iPeimbert et"all ([200I obtained Yp = 0.2477 ± 0.0029 
based on a sample of 5 H 11 regions, iPorter et al.l (|2005l ) 
He I emissiv ities, and adopting ii o n-zero temperature 
fluctuations. iFukugita fc Kawasakil (|2006D derived Yp = 
0.250 ± 0.004 for a sample of 31 H il regions studied by 
Izotov fc Thuaiij (|2004[ ). adopting iBeniamin et ail (|1999l 
He I emissivities. 
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It is generally believed that the accuracy of the de- 
termination of the primordial He abundance is limited 
presently, not so much by statistical uncertainties, but 
by our ability to account for systematic errors and bi- 
ases. There are many known effects we need to correct 
for to transform the observed He i line intensities into a 
He abundance. These effects are: (1) reddening, (2) un- 
derlying stellar absorption in the He i lines, (3) collisional 
excitation of the He i lines which make their intensities 
deviate from their recombination values, (4) fluorescence 
of the He i lines which also make their intensities deviate 
from their recombination values, (5) collisional and (6) 
fluorescent excitation of the H lines, (7) the temperature 
structure of the H 11 region and (8) its ionization struc- 
ture. All these corrections are at a level of a few percent 
except for effect (3) that can be much higher, exceeding 
10% in the case of the He i A5876 emission line in hot and 
dense H 11 regions. Most o f these effect s were analyzed 
and taken into account by llzotov et al.l (|2007| ) for their 
HeBCD sample. 

We present here a new determination of the primor- 
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dial He abundance. A new study is warranted because 
there are several recent new developments that allow 
more accurate estimates of some of the systematic effects 
mentioned above. Thus, iGonzalez Delgado et al.l ()2005l ) 
have calculated evolutionary stellar population synthesis 
models with high spectral resolution, so that equivalent 
widths of H I and He i absorption lines for different ages 
of single stellar populations and for a wi de range of metal- 
licities are now available. Furthermore. iLuridianal ()2009f ) 
has estimated the collision al enhan c emen t of hydrogen 
Balmer lines, while Luridia na et al.l (|2009f ) have consid- 
ered the fluorescent excitation of Balmer lines in gaseous 
nebulae (their case D). We incorporate all these new cal- 
culations in our new determination of the primordial He 
abundance. 

In 121 we briefly discuss the method used to derive He 
abundances in individual objects and the primordial He 
abundance from the total sample. In f|3l we derive the 
best value for Yp and the linear regression slope dY/dZ. 
The cosmological implications of our new results are dis- 
cussed in 21 We summarize our conclusions in Sj5l 

2. THE METHOD 

2.1. Linear regressions 

As in our previous work (see llzotov et al.ll2007l and 
references therein), we determine the primordial He 
mass fraction Yp by fitting the data points in the Y 
— 0/H plane with a linear r e gressi o n line of the form 
(iPeimbert fc Torres-PeimbertI Il974 119761 : iPagel et al.l 
[l99l 

y = K + ^7?^^(o/H), (1) 



d(0/H) 



where 



Y = 



%(1 - Z) 
l + 4y 
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is the He mass fraction, Z is the heavy element mass 
fraction, y = (y+ + y^+)x ICF(Re+ +}ie'^+) is the He 
abundance, y+ = Hc+/H+ and = He^^/H"'" are re- 
spectively the abundances of singly and doubly ionized 
He, and /CF(He+-|-He^+) is the ionization correction 
factor for He. 

The slope of the Y - 0/H linear regression can be 
written as: 
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where O and Z are respectively the mass fractions of 
oxygen and heavy elements. 

We also ta ke into account the depletion of oxygen on 
dust grains, llzotov et al.l ()2006!) demonstrated that the 
Ne/0 abundance ratio for low-metallicity BCDs is not 
constant, but increases with increasing oxygen abun- 
dance. This effect is small, with AlogNe/0 = 0.1 when 
the oxygen abundance changes from 12-|-logO/H — 7.0 
to 8.6. We attribute such a change to oxygen depletion, 
and correct the oxygen abundance for it, using the log 
Ne/0 versus oxyge n abundance regression line found by 
llzotov et al] (|2006D . and assuming that depletion is ab- 
sent in galaxies with 12-|-logO/H — 7.0. 

To derive the parameters of the Y versus 0/H lin- 
ear regressions, we use the maximum-likelihood method 



(jPress et al.lll99^ which takes into account the errors in 
both Y and 0/H for each object. 

2.2. A Monte Carlo algorithm for determining the best 

value of y"*" 

FoUowing llzotov et al.l ()2007D . we use the five strongest 
He I A3889, A4471, A5876, A6678 and A7065 emission 
lines to derive the electron number density Afe(He+) and 
the optical depth t(A3889). The He i A3889 and A7065 
lines play an important role because they are particularly 
sensitive to both quantities. Since the He i A3889 line is 
blended with the H8 A3889 line, we have subtracted the 
latter, assuming its intensity to be equal to 0.107 /(H/3) 
(jAUer 1984). 

The derived y"*" abundances depend also on a num- 
ber of other parameters: the fraction A/(Hq;)//(Hq;) 
of the Ha emission line flux due to collisional 
excitation, the electron number density iVe(He"''), 
the electron temperature re(He"'"), the equivalent 
widths EWafc«(A3889), EWa63(A4471), EWa63(A5876), 
EWahs(A6678) and EWats(A7065) of He i stellar absorp- 
tion lines, and the optical depth t(A3889) of the He i 
A3889 emission line. To determine the best weighted 
mean val ue of y"*" , we use the Monte Carlo procedure de- 
scribed in llzotov et al.l (|2007f ). randomly varying each of 
the above parameters within a specifled range. 

Additionally, in those cases when the nebular He ii 
A4686 emission line was detected, we have added to y"*" 
the abundance of d oubly ionized helium y^+ = He^^/H"*" 
(|Izotov et al.ll2007n . 

2.3. Parameter set for the He abundance determination 

For the determination of He abundance, we adopt He 
I emissivities from Port er et al.l (j2005) and take into ac- 
count the following systematic effects: 1) reddening; 2) 
the temperature structure of the H ii region, i.e. the 
temperature difference between Te(He+) and Te(0 ill); 
4) underlying stellar He i absorption; 4) collisional and 
fluorescent excitation of He i lines; 5) collisional and flu- 
orescent excitation of hydrogen lines; and 6) the ioniza- 
tion structure of the H ii reg i on. M ost of these effects 
were analyzed by llzotov et al.l ()2007D (see also references 
therein). 

We define the following s et of parameters : 

1. The r eddening law of iWhitfordI (1958) is adopted. 
llzotov et~a l. (2007) have shown that He abundances are 
not sensitive to the pa rticular redden i ng law adopted. 
For example, use of the iCardelli et al.l (|1989f ) reddening 
curve results in He and othe r element abund ances simi- 
lar to those obtained with the lWhitfordI ()1958l ) reddening 
law. The extinction coefficient C(H/3) is derived from the 
observed hydrogen Balmer decrement, after correcting 
the Ha, H/3, H7 and US line fluxes for the effects of col- 
lisional and fluorescent excitation. Finally, all emission 
lines are corrected for reddening, adopting the derived 
C(H/3). 

2. The electron temperature of the He+ zone is varied 
in the range Te(He+) = (0.95 ~ 1.0)xr^.(O ill). We have 
chosen this range following the work of iGuseva et al.l 
()2006f) and IGuseva et all ()2007D who have derived the 
electron temperature in the H"*" zone from the Balmer 
and Paschen discontinuities in the spectra of more than 
100 H II regions, and showed that Te(H"'") differs from 
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Te(0 III) by not more than 5%. We also assume that 
re(He+) = re(H+) because the H+ and He+ zones in 
our objects are nearly coincident. 

3. Oxygen abundances are calculated by considering 
two possible values of the electron temperature: 1) Tg = 
re(He+) and 2) = Te{0 ill). 

4. A^e(He+) and r(A3889) are varied respectively in the 
ranges 10 - 450 cm"'^ and 0-5, typical for extragalactic 
H II regions. 

5. The fraction of Ha emission due to coUisional excita- 
tion is varied in the range 0% - 5%, in ac c ordan ce with 
iStasihska fc Izotovl ()2001[ ) and iLuridianal H2OQ90 . The 
fraction of H/3, H7 and HJ emission due to collisional 
excitation is adopted t o be 60% that o f the Ha emis- 
sion, in acc ordance with ILuridianal (|2009f) . We note that 
llzotov eF^ l. (2007) underestimated that fraction for H/3, 
adopting a value of only 1/3, and neglected altogether 
to correct the H7 and MS emission lines for collisional 
excit ation. 

6. ILuridiana et al.l (|2009f ) have shown that the frac- 
tion of H/3 emission due to fluorescent excitation by the 
far-UV non-ionizing stellar continuum could be as high 
as 2%, and somewhat lower for the Ha emission (their 
case D). We have adopted the conservative value of 1% 
for the fraction of Ha, H/3, H7 and H5 emission due to 
fluorescent excitation, since a similar effect could affect 
the He i emission lines and partly compensate the effect 
for the Balmer H lines (the He abundance is calculated 
relative to that of H). 

7. The equivalent width of the He i A4471 absorp- 
tion li ne is chosen to be EW abs(A4471) = 0.4A, fol- 
lowing llzotov et al.l ()2007D and iGonzalez Delgado et al.l 
(|2005f ). The equivalent widths of the other ab- 
sorption lines are fixed according to the ratios 
EWafcs(A3889) / EWah.(A4471) = 1.0, EWa6.(A5876) 
/ EW„fc,(A4471) = 0.8, EW,6,(A6678) / EW,6,(A4471) 
= 0.4 and EWa6s(A7065) / EWa{,s(A4471) = 0.4. The 
EWa6s(A5876) / EWab.(A4471) and EWafas(A6678) / 
EWa6s(A4471) ratios were set equal to the values pre- 
dicte d for these ratios by a Starburst99 (jLeitherer et al.l 
Il999| ) instantaneous burst model with an age 3-4 Myr 
and a heavy element mass fraction Z = 0.001 - 0.004, 0.8 
and 0.4 respectively. These values are significantly higher 
th an the c orresponding ratios of 0.3 and 0.1 adopted 
by llzotov et al. (2007). We note that the value cho- 
sen for the EWa6s(A5876) / EWfl^^{A4£nj_raiioJs_aI^ 
consis tent with the one given bv IGonzalez Delgado et ahl 

Since the output high-resolution spectra in Star- 
burst99 are calculated only for wavelengths < 7000A, 
we do not have a prediction for the EWa&s(A7065) / 
EWa6s(A4471) ratio. We set it to be equal to 0.4, the 
value of the EWa6s(A6678) / EWa&s(A4471) ratio. 

8. The He ionization correction fact or 
/CF(He+-hHe++) is adopted from llzotov eTall (pOOTi ). 

3. THE PRIMORDIAL HE MASS FRACTION Yp AND THE 
SLOPE DY/DZ 

Two Y - 0/H linear r egres sions for the HeBCD galaxy 
sample of llzotov et al.l (|2007l ). with the above set of pa- 
rameters, are shown in Fig. [TJ The two regression lines 
differ in the way oxygen abundances have been calcu- 
lated. For the first regression line (Fig. [T^), oxygen 
abundances have been derived by setting the tempera- 
ture of the 0++ zone equal to T'e(He+), while for the 



second (Fig. [T|d), they have been derived by adopt- 
ing the temperature Te(0 ill) derived from the [O ill] 
A4363/(A4959-t-A5007) line flux ratio. 

The primordial values obtained from the two regres- 
sions in Fig. m = 0.2565 ± 0.0010 and Yp = 0.2560 ± 
0.0011, are very similar but are significantl y higher than 
the value Yp = 0.2516 ± 0.0011 obtained bv llzotov et~all 
(I2OO7I) for the same galaxy sample. The 2% difference 
is due to the inclusion of the correction for fluorescent 
excitation of H lines, the correction for a larger correc- 
tion for collisional excitation to the H/3 flux and larger 
adopted equivalent widths of the stellar He i 5876, 6678 
and 7065 absorption lines. We adopt the value of Yp from 
Fig. [1^, where both 0/H and Y are calculated with the 
same temperature = Te(He"'"). 

We have varied the ranges of some parameters to study 
how the value of Yp is affected by these variations. We 
have found that varying the fraction of fluorescent exci- 
tation of the hydrogen lines between 0% and 2%, and/or 
setting Te(He+) = TeiO ill) or changing Te(He+) in 
the range (0.9 - 1.0) x Te{0 III) (instead of making it 
change between 0.95 and 1.0 x Te(0 in)), result in a 
change of Yp between 0.254 and 0.258. Additionally, 
adding a systematic error of 1% caused by unc ertain- 
ties in the He i emissivities (jPorter et al.l l2009f) , gives 
Yp^ 0.2565 ± 0.0010(stat.) ± 0.0050(syst.), where "stat" 
and "syst" refer to statistical and systematic errors, re- 
spectively. Thus, the value of Yp derived in this paper, 
is 3.3% greater than the value of 0.2482 obtained from 
the 3yr WMAP data, assuming SBBN (Sperg cl et al.l 
120071 ). Howe ver, it is consist e nt wit h the Yp = 0.25^°;^° 
obtained by 'I chikawa et al.l (2008') from the available 
WMAP, ACBAR, CBI, and BOOMERANG data [actu- 
ally, the peak value in their one-dimensional marginalized 
distribution of Yp (their Fig. 3) is equal to 0.254]. 

Using Eq. |3l we derive from the Y - 0/H linear re- 
gression (Fig. [TJi) the slopes dF/dO = 2.46±0.45(stat.) 
and dY/dZ = 1.62±0.29(stat.). These slopes are shal- 
lo wer than the o nes o f 4.33±0.75 and 2.85±0.49 derived 
bv llzotov et"all (|2007D . 

4. DEVIATIONS FROM SBBN 

We now use our derived value of the primordial He 
abundance along with the observed primordial abun- 
dances of other light elements to check the consistency 
of SBBN. Deviations from the standard rate of Hubble 
expansion in the early Universe can be caused by an ex- 
tra contribution to the total energy density, for example 
by additional flavors of neutrinos. The total number of 
different species of weakly interacting light relativistic 
particles can be conveniently be parameterized by N^,, 
the "effective number of light neutrino species" . 

To perform the study, we use the statistical tech - 
niq ue, with the code described bv lFiorentini et all ()1998D 
and ILisi et all (1999.) . This code allows to analyze the 
constraints that the measured He, D and ^Li abun- 
dances put on 77 and N^. For the primordial D abun- 
dance, we use the value obtained bv lPettini et all ([2008). 
As for ^Li, its value derived from ob servations of low- 
metallicity halo stars (|Asplund et al.ll2006 ') is ~ 5 times 
lower than the one o btained from the WMAP analysis 
(jDunklev et all l2009f ) . Because mechanisms that may 
lead to a reduction of the ^Li primordial abundance, such 
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as diffusion or rotationally induced mixing, are not well 
understood and we do not know how to correct for them, 
we have adopted the value of the primordial abundance 
of ^Li abundance as de rived from the Syr WMAP data 
of lDunklev et alJ (|2009') . The predicted primordial abun- 
dances of light elements depend on the adopted neutron 
life-time t^. We have considered two values, th e old one, 
r„ = 885.4 ± 0.9 s (lArzumanov et all I2000D. a nd the 
new one, t„ = 878.5 ± 0.8 s (jSerebrov et al.ll200l 12001 . 
With the old value of Tn and our best value of the pri- 
mordial He abundance, Yp = 0.2565± 0.0010(stat.) ± 
0.0050(syst.), the minimum xLin (= 0.640524) is ob- 
tained when riio = 6.47 and N^, = 3.68. The value of 
7710 is in agreemen t with tjiq = 6 . 23ib0 .17 derived from 
the WMAP data (IPunklev et"al] I2009D . If instead the 
new value of is adopted with the same value of Yp, 
then the minimum Xmin i~ 0.619816) is obtained when 
Vio — 6.51 and Ni, = 3.80. We note that 7710 and N,y only 
slightly depend on the value of the ^Li abundance. They 
are decreased by 3% a nd 2% , respectively, if the observed 
^Li abundance by 1 Asplun d et al. (2006) is adopted. 

The joint fit of 77 and N^, is shown in Figures [2^ and 
[2}d for the two values of Tn. The la (x^ - Xmin — 1-0) 
and 2a (x^ - Xmin — 2-71) deviations are shown respec- 
tively by the thin and thick solid lines. We find the 
equivalent number of light neutrino species to be in the 
range = 3.68tE;-|o {2a) (Fig. Hi) in the first case, 
and = 3.80l;o|E] {2a) (Fig. Eb) in the second case. 
Both of these values are only marginally consistent (at 
the 2a level) with the experimental value of 2.993±0.011 
(|Caso et al.l I1998D shown by the dashed line, implying 
deviations from SBBN. We note that, although both val- 
ues are consistent with N^, = 4.4 ± 1.5 derive d from 
the a nalysis of Syr WMAP observations (Komats u et al.l 
12009^) , the primordial helium abundance sets tighter con- 
straints on the effective number of neutrino species than 
the CMB data, the error bars of N^, being approximately 



half as large in the first case as compared to the latter 
case. 

5. SUMMARY AND CONCLUSIONS 

We present here a new determination of the primordial 
helium mass fraction Yp by linear regressions of a sam- 
ple of 93 spectra of 86 low-metallicity extragalactic H 11 
regions. 

In this new determination of Yp , we have taken into ac- 
count the latest developments concerning several known 
systematic effects. We have used Monte Carlo methods 
to solve simultaneously for the effects of collisional and 
fluorescent enhancements of He i recombination lines, of 
collisional and fluorescent excitation of hydrogen emis- 
sion lines, of underlying stellar He i absorption, of possi- 
ble temperature differences between the He"*" and [O ill] 
zones, and of the ionization correction factor /CF(He"'" 
+ He2+). 

We have obtained the following results: 

1. Our best value is Yp = 0.2565 ± 0.0010(stat.) ± 
0.0050(syst.), or 3.3% larger than the value derived from 
the microwave background radiation fluctuation mea- 
surements assuming SBBN. In order to bring this high 
value of Yp into agreement with the deuterium and 5yr 
WMAP measurements, an equivalent number of neutrino 
flavors in the range 3.68 - 3.80, depending on the life- 
time of the neutron, is required. This is higher than the 
canonical value of 3 and implies the existence of devia- 
tions from SBBN. 

2. The dY/dZ slope derived from the Y -O/R linear 
regression is equal to 1.62 ± 0.29(stat.), shallowe r than 
the previous determination bv llzotov et al.l (|2007f ). 
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Fig. 1. — Linear regressions of the helium mass f ractio n Y vs. oxygen abundance for H II regions in the HeBCD sample. The Ys are 
derived with the He I emissivities from lPorter et aH 120051) . The electron temperature Te(He+) is varied in the range (0.95 - l)xTe(0 III). 
The oxygen abundance is derived adopting an electron temperature equal to Te(He+) in a) and to Te{0 III) in b). 
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Fig. 2. — a) Joint fits to the bary on-to-photon number ra tio, ?7 io, and th e equiv alent number of light neutrino species A'^^, using a 
analysis with the code developed by ll^'iorentini et al] 119981) . and |Lisi et a Tl p999l) . The value of the primordial He abundance has been 
set to Yp = 0.2565 (t his paper), that of (D/H)p is taken from[Pettini et al. (2008) and that of (^Li/H)p from 5yr WMAP measurements 
IDunklev et al.l 120091 ). A neutron lifetime Tn = 885.4 ± 0.9s from Ar zunianov et al. (200 0.) has bee n adopted. Thin and thick solid lines 
represent respectively la and 2(t deviations. The experimental valu e Ni, = 2 .993 (Caso et al.lll998l) is shown by the dashed line, b) The 
same as m (a), but with a neutron lifetime Tn = 878.5 ± 0.8s (.Serebrov et al.l l2005. 2008). 
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